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Abstract. Some recent developments concerning the role of strange quark matter 
for astrophysical systems and the QCD phase transition in the early universe are 
addressed. Causality constraints of the soft nuclear equation of state as extracted 
from subthreshold kaon production in heavy-ion collisions are used to derive an upper 
mass limit for compact stars. The interplay between the viscosity of strange quark 
matter and the gravitational wave emission from rotation-powered pulsars are outlined. 
The flux of strange quark matter nuggets in cosmic rays is put in perspective with a 
detailed numerical investigation of the merger of two strange stars. Finally, we discuss 
a novel scenario for the QCD phase transition in the early universe, which allows for 
a small inflationary period due to a pronounced first order phase transition at large 
baryochcmical potential. 



1. Introduction 

There is no attempt in this contribution to discuss all the recent advances in the very 
active research area of studying the impacts of strange quark matter in astrophysics 
and cosmology. We refer to the other contributions of this conference for topics not 
covered in this contribution to the proceedings. Signals for the QCD phase transition 
to strange quark matter in core-collapse supernovae are addressed in the contribution 
by Irina Sagert et al. (see also [I]). For the appearance of a new strange phase in the 
evolution of proto-neutron stars we refer to the contribution by Giuseppe Pagliara et 
al. (see also j2]). The issues of how to nucleate strange quark matter in astrophysical 
systems has been studied in detail in the contribution of Bruno Mintz et al. (see also 
[3]). Strange quark matter, colour superconductivity and the relation to quarkyonic 
matter are discussed in the contribution by David Blaschke et al. (see also [1]). 

In the following we give attention to a few other recent developments in the field. 
First we reconsider the maximum mass constraint for neutron stars using causality 
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arguments for the nuclear equation of state. Constraints on the nuclear equation of 
state using heavy-ion data, here subthreshold production of kaons, can be adopted to 
give a new maximum mass limit of about 2.7 solar masses. Then we are studying the 
gravitational wave emission from rotating neutron stars and the relation to the viscosity 
of strange quark matter. The recent simulation of the merger of two strange stars, the 
release of strange quark matter nuggets to the interstellar medium and the possible 
impact for the flux of strangelets in cosmic rays are another issue being highlighted. 
Finally, we are concerned with the cosmological QCD phase transition which happens 
around 10 -5 seconds after the big bang. If one allows for an inflationary period at a 
strong first order phase transition, the early universe can actually pass through the QCD 
phase diagram at large baryochemical potential in contrast to the standard model of 
cosmology. 

2. Maximum mass of neutron stars: constraints from heavy-ion data 

Neutron stars are produced in core-collapse supernova explosions and are compact, 
massive objects with radii of 10-15 km and masses of 1-2 M Q . Pulsars are rotation- 
powered neutron stars and have been observed in binary systems, i.e. with a companion 
star, with a white dwarf or even with another neutron star. To create such a binary 
pulsar one starts with two ordinary stars, one with at least 8M . The heavier star 
explodes first in a supernova and leaves a neutron star as a remnant. The neutron star 
is spun up by accreting matter from the companion star. Then also the companion star 
might become a white dwarf or another neutron star depending on its initial mass (see 
[5] for a review). 

For binary pulsars masses can be pretty nicely determined. At present more 
than 1800 pulsars are known with 140 binary pulsars. From those binary pulsars, 
the pulsar PSR J1903+0327 is particular interesting [61 |7]. By measuring the post- 
Keplerian parameters from pulsar timing, here the Shapiro delay parameters r and s 
alone, constrain the mass to M = (1.67 ± O.11)M . Combined with the periastron 
advance Cj the mass is constraint to be M = (1.67 ± O.O1)M . This mass measurement 
would put some severe constraints on soft nuclear equations of state which give a smaller 
maximum mass [8]. 

Kaons are produced by associated production: NN— > NAK and NN— >-NNKK in 
elementary collisions. In the medium as generated in central heavy-ion collisions new 
processes emerge through rescattering as 7rN — > AK, 7rA — > NK so that kaons can be even 
produced below the (elementary) threshold bombarding energy. In heavy-ion collisions 
of 1-2 AGeV nuclear matter can compressed up to 3no- As kaons have a long mean- free 
path, they can escape the high density region serving as a messenger of the properties 
of high-density nuclear matter (for a review see e.g. [9]). One defines the double ratio 
as the multiplicity of produced kaons per mass number for Au+Au collisions relative to 
the one for C+C collisions. This observable turns out to be rather insensitive to input 
parameters but strongly dependent on the nuclear equation of state (to be more precise 
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the in-medium potential of nucleons which can be related to the nuclear equation of 
state). Using a simple Skyrme parametrization, only simulations with a compression 
modulus of Kn ~ 200 MeV [THl E] can describe the energy dependence of the kaon 
production data as measured by the KaoS collaboration [T2J [13]. Hence, the nuclear 
equation of state would be soft around 2-3 uq. 

Causality arguments can then be adopted to constrain the maximum possible mass 
of neutron stars. Above a fiducial density as determined from the data analysis of the 
KaoS data the stiff est possible equation of state is taken which gives the maximum 
pressure for a given energy density allowed. Causality demands then that at most 
p = e— e c above the fiducial density e/ (see [T^lfTB]). There is a general scaling relation for 
the maximum possible mass of compact stars which is given by M max = 4.2M (eo/e/) 1 ^ 2 
where €q is the energy density of normal nuclear matter (the numerical prefactor depends 
slightly on the low-density nuclear equation of state, see [16]). For a fiducial density 
of €f = 2e as provided by the heavy- ion data the upper mass limit is about 2.7 M Q 

pang. 

3. Gravitational wave emission from rotating neutron stars: the viscosity of 
strange quark matter 

There exists a special class of binary neutron stars, so called x-ray burster. In these 
astrophysical systems, the neutron star is accreting matter from an ordinary star or 
a white dwarf. A low-mass x-ray binary (LMXB) is a neutron star with a low mass 
companion star or a white dwarf. The accreting matter falls onto the neutron star 
surface and thereby releases x-rays in a burst. The neutron star surface is heated up to 
typical temperatures of T = 10 — 100 keV. 

For rotating neutron stars, there exists an instability, the r-mode instability. 
Oscillations brings the neutron star matter out of /3-equilibrium which is restored by 
weak processes affecting the bulk and shear viscosity [19] . The viscosity of nuclear matter 
is quite low in certain regions of frequency and temperature so that the oscillations will 
not be damped, gravitational waves are emitted and the star's rotation slows down. The 
contributions from quarks increase the stability window in temperature and frequencies 
substantially. In color-superconducting phases the viscosity of quark matter can be 
drastically different [20, EH [221 E3], EI] . In the 2SC phase, with the pairing of two quark 
flavours and two colours, there are still unpaired quarks left and one obtains similar 
results compared to the free case [25, 26J. Both cases are compatible with the data on 
low-mass x-ray binaries. On the contrary, quarks of all color and flavor are paired in the 
CFL phase. The reaction rates are suppressed, the instability window is considerably 
broader which could leave an observable impact on the evolution of young pulsars and 
accreting neutron stars. Recently, it was also found that there is a new process which 
takes away the energy from the r-mode oscillations by the emission of neutrinos. It 
acts effectively as a viscosity which would help in damping the r-mode instability by a 
substantial factor [27]. 
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4. Strange star mergers: the flux of strangelets in cosmic rays 

According to the Bodmer-Witten hypothesis absolutely stable strange quark matter 
would be more bound than ordinary nuclear matter. If one allows for absolutely stable 
strange quark matter, then strange stars will exist which are compact stars consisting of 
absolutely stable strange quark matter being bound by strong interactions (selfbound) 
and not by gravity as ordinary neutron stars. Strange stars can coexist with neutron 
stars, but in principle a neutron star would collapse to strange star in a spectacular 
astrophysical event. Recently, the collisions of two strange stars have been simulated 
in 3D relativistic smoothed particle hydrodynamics with approximate treatment of 
effects from General Relativity [28]. Typical timescales involved are milliseconds, so 
that matter is in /3-equilibrium. The simulations have been performed for two bag 
constants of B = 60 and 80 MeV/fm~ 3 , in both cases one gets absolutely stable strange 
quark matter, and for different initial masses of the two strange stars. In principle, two 
different scenarios have been found: either there is an immediate collapse to a black hole 
or a hypermassive object stabilized by differential rotation is formed. The implications 
of these results are quite striking: the ejected mass of strange quark matter (strangelets) 
released to the interstellar medium are estimated to be (1.4 — 2.8) ■ 10 _4 M Q for the lower 
bag constant and zero for the larger one. As strange star mergers are considered to be 
the prime source for strangelets in cosmic rays, the latter result opens the possibility 
that no strangelets are seen in cosmic rays although strange stars exist. The Alpha 
Magnetic Spectrometer AMS will measure the strangelet flux in cosmic rays (planned 
to be installed at the International Space Station ISS in 2010). Interestingly, the 
gravitational wave pattern of compact star mergers can be used to distinguish between 
neutron star and strange star mergers [29] . 

5. Early universe: phase transition from strange quark matter 

Let us briefly outline the history of the early universe. In a radiation dominated universe 
the temperature scales inversely with the scale parameter a. Big-bang nucleosynthesis 
happens at times of t = Is to 3 minutes which corresponds to temperatures of T = 1 to 

0. 1 MeV. The QCD phase transition occurs at a temperature of about T w 170 MeV, 

1. e. at a time of t ~ 10~ 5 s after the big bang. The early universe passes through the 
electroweak phase transition at t 10~ 10 s (T « 100 GeV). In the standard cosmology, 
one knows from the measurement of the microwave background radiation with WMAP 
and from the observed light element abundance from the big bang nucleosynthesis that 
the ratio of the number of baryons to the number of photons (which is related to the 
total entropy of the universe) is tiny: 

Note that the baryon number per entropy is a conserved quantity throughout the 
evolution of the universe unless there is a non-equilibrium process as a first order phase 
transition. From the above considerations one concludes that the early universe evolves 
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Figure 1. The evolution through the QCD phase diagram for the little inflation 
scenario. The universe is trapped in a metastable vacuum at the first order phase 
transition line, supercools and is reheated back close to the critical temperature at 
tiny values of the baryochemical potential. 



along fx ~ which implies that the early universes crosses the QCD phase transition 
in a region of the QCD phase diagram where one expects a crossover transition as 
seen in recent lattice gauge calculations. A crossover transition implies that there is 
nothing spectacular happening, as one expects no cosmological signals from such a kind 
of transition. 

The Friedmann equation for radiation dominated universe reads 

tt2 8tvG _ T 4 
H^— p ~ gi T)- 

where g(T) is effective number of relativistic degrees of freedom at a given temperature 
T and M p is the Planck mass. The Hubble time is defined by the Hubble parameter 
and is related to the true time by t = 3tn for a radiation dominated universe. By using 
the Friedmann equation one finds that tn = l/H ~ g~ l l 2 Mp /T 2 so that there is a 
simple relation between the time after big bang and the temperature of the radiation 
dominated early universe 

t /!MeV\ 2 



') 



1 sec V T 

We see that if the number of degrees of freedom changes smoothly with the temperature, 
as it would be the case for a crossover phase transition, the temperature changes 
continuously with time. 

The interesting question is what would happen if the early universe passes through 
a first order phase transition? We argue in the following that this is a viable scenario 
for the QCD phase transition and not in contradiction to our present knowledge of 
QCD. Furthermore, we demonstrate that it would leave an observable imprint on the 
gravitational wave background [30] . 

A first order phase transition in general allows for a false metastable vacuum state 
where the universe could be trapped for some time. For a constant vacuum energy, 
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one arrives at the de Sitter solution of the Friedmann equations which results in an 
(additional small) inflationary period: 

H = a I a ~ M~ l p\l 2 — H v — const. — > a ~ exp(H v ■ t) 

Remarkably, just a few e-folds are enough (standard inflation needs iV ~ 50) for 
our purposes. The baryon to photon ratios before and after the inflationary period 
should scale with the ratio of the scale parameters cubed as entropy is produced during 
reheating while baryon number is conserved so that 




The final ratio should be 10 -9 as observed so that we need just a boost of N = 
In (ay/aj) ~ ln(10 3 ) ~ 7, i.e. seven e-folds, to get an initial ratio of (fi/T) { ~ 0(1). 
Note that for standard inflation one needs about N = 50 e-folds, so that the little 
inflation at the QCD phase transition can not replace it but could be present as an 
additional small inflationary period in the early universe. 

So our scenario is as follows: the early universe is initially at large baryochemical 
potentials ji/T ~ 1. Such a large value is possible for e.g. Affleck-Dine baryogenesis 
which involves scalar supersymmetric fields carrying baryon number. The early universe 
reaches the first order phase transition line of QCD at high baryochemical potentials 
and is trapped in the false vacuum. Now the inflationary period starts with supercooling 
and dilution with fi/T = const.. The decay to the true vacuum state will release latent 
heat so that the universe is reheated to T ~ T c . Afterwards, the final baryon to photon 
ratio is given by ji/T ~ 10~ 9 due to the entropy produced during the transition. Finally, 
the universe evolves along the standard cosmological path to big-bang nucleosynthesis 
and so on. The path through the QCD phase diagram is depicted in Fig. [TJ 

However, there are observable differences to the standard model of cosmology. The 
first order phase transition produces perturbations, scalar and tensor ones, by colliding 
bubbles. The tensor modes correspond to gravitational waves with a frequency scale 
given by the (redshifted) horizon scale at the transition point 

^peak ~ H ■ T 7i0 /T QCD ~ T QCD /M p • T 7j0 ~ 1(T 7 Hz 

The maximum amplitude of the gravitational waves is then given by the ratio of the 
scale factors today and at the QCD phase transition and is h ~ oqcd/^o ~ 10~ 12 . For 
superhorizon modes v < H the amplitude scales as h(u) oc v~ x l 2 which corresponds 
to white noise. The modes within the horizon v > H are determined by multi bubble 
collisions and could scale as h(u) oc v n where n is between —1 and —2 J3TJ [32]. The 
gravitational waves from a first-order QCD phase transition are observable with pulsar 
timing, as with the Parkes Pulsar Timing Array and in the future with the Square 
Kilometre Array SKA (see figure [2]). 

Other cosmological implications of a first-order QCD phase transition are that the 
cold dark matter density is diluted by the same factor as the baryon number density, 
i.e. by up to a factor 1CT 9 . To explain the present cold dark matter relic density within 
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Figure 2. The gravitational wave spectrum from a first order QCD phase transition. 
The Parkes Pulsar Timing Array PPTA already gives limits on the maximum 
amplitude of those gravitational waves. In the future the Square Kilometre Array 
SKA will push these limits down by a few orders of magnitude. For a flat spectrum at 
high frequencies, LISA would be also sensitive. 



the standard WIMP scenario, one needs a different WIMP annihilation cross section 
o"ann as the cold dark matter density scales as ^cdm ~ Cweak/oann or by a substantially 
larger WIMP mass. Both cases can be probed by the LHC! The large-scale structure 
power spectrum responsible for galaxy and galaxy cluster formation would be modified 
up to a mass scale of M ~ 10 9 M Q which is the typical mass scale of dwarf galaxies. 
Without QCD inflation, effects can be only up to the horizon mass ~ 1O~ 9 M at the 
phase transition point. Also, seeds of (extra)galactic magnetic fields can be generated 
by charged bubble collisions during the phase transition. The observed magnetic field 
today in our galaxy is B ~ IO -5 G, extragalactic ones are in the range of B ~ IO -7 G. 
One needs primordial seed fields of B — 10~ 30 . . . 10~ 10 G. As the electroweak phase 
transition is a crossover for large Higgs masses, a first order QCD phase transition 
in the early universe would provide then a possible explanation of those cosmological 
magnetic fields within the standard model again (see |33j for a detailed calculation). 
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